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SUMMARY 

An analytical and experimental investigation was made of the 
effect of twist on the vibrations of cantilever beams. The analyt- 
ical investigation was made by the use of Station Functions. Gen- 
eral equations are developed for the coupled bending - bending- 
torsion vibrations of a cantilever beam, and it is shown how these 
equations reduce in simpler cases. The use of tabulated Station 
Numbers makes tjie analytical method presented herein particularly 
simple to apply. 

An example is worked out in detail to illustrate the applica- 
tion of the method. Good agreement is obtained among the method 
presented, an exact theoretical solution developed herein, and 
experimental results. 

It is shown that for a beam with a ratio of bending stiffness 
in the two principal directions equal to 144, the effect of coupling 
due to twist is to raise the value of the first natural frequency 
by a negligibly small amount, to decrease steadily the second fre- 
quency, and to lower the third frequency considerably. 


INTRODUCTION 


The failure of turbine and compressor blades due to vibrations 
has led to an increased interest in the study of the vibrations of 
these blades and in the determination of the natural modes and fre- 
quencies. In such theoretical studies, the compressor or turbine 
blade is usually assumed to act as a cantilever beam. The calcula- 
tion of the uncoupled modes of arbitrarily shaped cantilever beams 
has been extensively investigated, but little work has as yet been 
done on calculating the coupled modes of such beams. If the geometry 
of the beam is such that coupling exists, the coupled modes are 
the, actual vibrational modes that must be calculated. 
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There are three types of coupling that generally need to he 
considered: (1) coupling between bending and torsion, (2) coupling 

between bending in two directions due to natural twist, and (3) 
coupling between bending in two directions and torsion. Bending- 
torsion coupling has been studied by several investigators (ref- 
erences 1 to 4), but little work has been done on coupling involv- 
ing twist. 

The earliest work on the vibrations of twisted beams involved 
propeller blades and the general practice was to neglect the effect 
of twist. Later investigations include the effect of twist by 
applying the methods of static influence coefficients (references 5 
and 8) or integral equations (references 6 and 7). 

In the use of integral equations, great difficulty is encountered 
in the calculation of the higher modes of vibration because exact 
orthogonality conditions are required to prevent convergence back 
to the lowest mode. The influence coefficient method reduces the 
problem to one having a finite number of degrees of freedom. _ The 
accuracy of the higher modes, however, is poor. In using such a 
finite number of degrees of freedom, only the first third of the 
modes and first half of the frequencies obtained are within the 
usual engineering accuracy. 

The present investigation was made at the NACA Lewis laboratory 
in order to present a straightforward accurate method for determin- 
ing the coupled modes and frequencies of nonuniform twisted canti- 
lever beams and to determine the effect of twist on the vibrational 
frequencies of such beams. The method is based on the use of Sta- 
tion Functions as developed for uncoupled and coupled bending- 
torsion vibrations in reference 4 and first discussed in refer- 
ence 9. Incorporated in the method are the advantages of the 
continuous-function deflections of the Rayle igh-Ritz and Stodola 
methods together - with the advantages of the finite number of 
degrees of freedom of the influence- coefficient method. Use is 
made of the Station Numbers derived in reference 4 and tabulated 
herein for convenience. An example is given and comparisons are 
made with exact theoretical and experimental values. 


ANALYSIS 

The present analysis of the vibrations of twisted cantilever 
beams is carried out in two ways: (1) an approximate solution 

using Station Functions and (2) an exact solution of the differen- 
tial equations of equilibrium of the system. 
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The approximate solution using Station Functions for the gen- 
eral case of coupled tending - bending-torsion vibrations is pre- 
sented herein. It is shown how the general equations reduce to 
the special cases of bending-bending, bending- torsion, bending, and 
torsion vibrations. 

All symbols used in this section are listed and defined in 
appendix A. The exact solution and the derivation of the Station 
Function equations presented are discussed in appendixes B and C. 

In the present analysis, matrix notation is used entirely. 

The use of matrix algebra reduces the amount of labor considerably. 
The only matrix operations needed are multiplication and addition. 
For those unfamiliar with the multiplication and addition of 
matrices, these simple processes are described in appendix D. The 
use of matrices, however, is not necessary. The scalar equations 
given in appendix B can be used for all the computations. 

Bending - bending- tors ion vibrations. - It is shown in appen- 
dix B that the deflections of the reference stations for a beam ' 
divided into n intervals of length 8, as shown in figure 1, are 
given by the following matric equation: 



o^hiqB 4 
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( 1 ) 


where the elements of these matrices are themselves matrices given 
as follows: 
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D 0 5 T 1 (KL + EM) 

D s T (K L’ + H M’) 

0y l y y 

D e x H T i ( K x L ’ + *X M ’> 

D = T 2 F P’ - T x F Q' + T 3 Gy N* + T 4 G y M' 

( 2 ) 


D, 


D yx “ D xy 


D. 


ye 


D. 


X0 


H T 2 F x P' - F x Q* + T 3 G x N’ + T 4 G x M' 

h T 2 F xy p- - T x F xy Q- + T 3 G xy N’ + T 4 G xy M* 

£ t 2 F ye p - T i F ye * + T s G ye N + t 4 Gye M 

£ T 2 F x0 P - T 1 F xe Q + T 3 G x0 W + T 4 G xe m 


The matrices L, M, H, P, Q, L', M', N’, P‘, and Q« are 
matrices of Station Numbers. These matrices are functions only 
of n, the number of stations, and have been tabulated in tables I 
to VIII for n equal 1 to 8. 


The matrices T,, Tg, T 3 , and T 4 are standard matrices listed 

in table IX. The first n rows and columns of these matrices are 
used, where n is the number of stations taken along the beam. 

The H, K, F, and G matrices involve the physical properties 
of the beam. The element cf the i th row and j th column of each' 
of these matrices is given as follows: 


Let 
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G y0i.j ( f yi S yj 

G x0ij = ^ f xi S xj 
The matrices f , Y, and X 


f xyi S xj^ 7 ij 

" f xyi S yj ) 7 ij 

are column matrices given fcy 
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The corresponding scalar equations, equations (B5) and (B6), 
c^n he found in appendix B. 

The labor involved in the calculation of the D matrices of 
equation (2) can be considerably reduced if the equations are 
expressed in the following partitioned form: 
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With 


becomes 


1 EI xO 
“J 2 m Q 8 4 


= X. , the characteristic equation for the system 
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or 

|d -xi | = o 

where D is the matrix indicated, the dynamical matrix of the sys- 
tem, and I is the identity matrix. The corresponding scalar 
equations (B7) can be found in appendix B. 

If the physical properties of the beam under consideration are 
known for each of the n intervals, the F, G, H, and K matrices 
denoting these physical properties can be written. The D matrices 
can be calculated by straightforward matrix computation using equa- 
tions (2) and (2a) and tables I to IX. Any method available can 
then be used to. solve equation (3) for the 3n values of 

Special cases . - For the cases where certain couplings between 
modes are negligible, equation (3) reduces to the following special 
cases: 
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3_ EIxO \ 

(1) Bending-bending vibrations. With -g j = \ 

co m^S 


v w 


yx 


D xy 

V^ 1 


= 0 


(2) Bending- tors ion vibrations. With 


1 EI xO _ 


uj 2 m 0 5 4 


= X 


(3a) 


Td-Ki e Td 


D. 


ye 


ey 

V XI 
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(3b) 


1 EE xO 

(3) Uncoupled bending vibrations. With — 

or n^S 4 


= X. 


|Dy- XI =0 


1 c o 

(4) Uncoupled torsion vibrations. With — r r 

“ V 


= X. 


(3c) 


D 0 - XI I = 0 


(3d) 


The determinantal equations (3) can be solved for the modes 
and the frequencies of the type vibration under consideration. 


EXAMPLES 

In applying the previously discussed method, the elements of 
the dy nami cal matrix must be determined for a given beam. These 
quantities depend upon the physical properties of the beam and upon 
the number of stations chosen. If the physical properties of the 
beam are known, the dynamical matrix can be calculated directly 
from- equations (2) or (2a). The matrices T-^, Tg, T^, T^_, L, M, N, 

P Q L' M' W’,P', and Q' appearing in equations (2) depend 
only upon the number of stations n used and can be read directly 
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from tables I to IX. Then, one of equations (3) can "be solved for 
the modes and the frequencies of vibration by any method desired. 

The following detailed example illustrates the method for coupled 
bending - bending- torsion vibrations: 

Consider a cantilever beam with I^/l|, the ratio of bending 

stiffness in the £ -direction to bending stiffness in the 
q-direction, equal to 100, and having linear twist, with total 

i Iq Sly 

angle of twist equal to 1 radian. Let T = — = 0.03, 

6 2 C 0 m 0 

let r^/rj: equal 10, and let € , the coefficient of coupling 

between torsion and bending -in the q-direction, r yQ 2 / r go 2 > 0< l ua - 1 - 0.2 . 

Let the number of stations chosen be two (n = 2). The calcu- 
lation of the elements of the F, G, H, and K matrices is shown 
in table X. The multiplication of these matrices by the appropri- 
ate L, M, N, P, Q, L', M», N», and Q' matrices from table II 
is also shown, and the determinantal equation is calculated. The 
roots of this 6 X 6 determinant can be computed in any manner 
desired. 


APPLICATIONS AND RESULTS 


The method was applied to determine the effect of twist upon 
the natural frequencies of a cantilever beam. Calculations were 
performed upon a beam such as that shown in figure 2. The 
ratio I^/l ^ was chosen as 144. Calculations were made for several 

values of total twist, ranging from 0 radians to 1.0 radian. The 
first, second, and third frequencies of vibration were obtained in 
each case. 


The procedure for obtaining the exact theoretical values is 
derived in appendix C. A determination of the first natural fre- 
quency was made for I^/l^ equal to 144, and having Cp^, the total 

angle of twist, equal to 0.70 radian. 


Several beams having length equal to 6^ inches and 

u 


yq - 144 

were machined, mounted, and vibrated to determine their natural 
frequencies. These frequencies were determined to within ±1 per- 
cent. Several values of total twist were taken. 
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The results of the calculations performed are plotted in fig- 
ure 3 for the first, second, and third modes of vibration. The . 
effect of coupling due to twist on the first mode for this beam is 
seen to raise the natural frequency to a value slightly higher than 
it would have if there were no twist; for a total angle of twist 
of 0.75 radian, the percentage increase is 0.5 percent. This dif- 
ference is less than the experimental error. 

Coupling due to twist has a more pronounced effect upon the 
higher modes of vibration. For the second mode, the natural fre- 
quency decreases to approximately 82 percent of its untwisted value 
for a total angle of twist equal to 0.75 radian and decreases to 
approximately 72 percent of its untwisted value for a total angle 
of twist equal to 1.0 radian. 

The natural frequency of the third mode also decreases. For 
a total angle of twist of 0.75 radian, the frequency is approxi- 
mately 80 percent of the untwisted value. For this mode, the per- 
centage decrease first becomes larger and then smaller as the total 
angle of twist is increased through the range investigated. 

The experimental results are also shown in figure 3. The 
beams tested had total angles of twist equal to 0.11, 0.39, 0.6, 
and 0.76 radian. The first mode changes very little. For the 
second mode, the percentages of the untwisted value for the fre- 
quency are 101, 93, 88, and 75 percent. The third- mode frequen- 
cies could not be obtained with any degree of certainty. 


In addition, experimental data taken from reference 10 have 
been interpolated for a value of X^/l^ = 144 and are plotted in 

figure 3. Agreement between the theoretical and experimental 
results, in general, is very good. 


CONCLUDING REMARKS 

Calculations based on the method presented tc determine the 
effect of coupling due to twist on the natural frequencies of. 
twisted cantilever beams were performed. For a beam with ratio of 
bending stiffness in the two principal directions equal to 144, 
the first natural frequency is raised by a negligibly small amount; 
for the second mode, the natural frequency decreases greatly as the 
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twist is increased; the effect of coupling due to twist for the 
third mode is also to lower its frequency considerably. 

A comparison among the method presented, an exact theoretical 
solution, and experimental results was made for a specific case; 
good agreement resulted. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, October 17, 1950. 
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APPENDIX A 


SYMBOLS * 

The following symbols are used in the report: 

torsional stiffness of beam at root 

ratio of average torsional stiffness of i 
val to torsional stiffness at root 


th 


inter- 


E 

fjUb SjU) 

f f f 

V y’ xy 
f xi> f yi> f xyi 


0 


7 * 


h 


’ \ 


x0 

i, j, k, r 

l 

m 


Young's modulus 

Station Function for torsion and bending, respec- 
tively, associated with j^h station 

I x /El£l v Iy/EI^, and I^/EI^, respectively 

average value of Ljj.I^.q/ 1^1^ , Iyl^/ l£ 1^ , and 
I X yI x0 /l^I T) over i th interval 

identity matrix 

mass moment of inertia per unit length of beam at 
root 

ratio of average mass moment of inertia per unit 
length of beam of i" t * 1 interval to mass moment 
of inertia per unit length at root 

moment of inertia of cross section of beam about 
x- , y-, £- , and p-axes, respectively, function 
of z 

moment of inertia of cross section of beam about 
x-axis at root 


summation indices 
length of beam 

mass per unit length of beam, function of z 
mass per unit length of beam at root 


mo 
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ratio of average mass per unit length of 

i^ h interval to mass per unit length at root 

n number of intervals 

r y , r x , r^ , r, absolute magnitude of projection of distance from 
5 elastic axis to center of gravity on perpendicu- 
lar to y, x > Ti, and { bending directions, 
respectively 

r y Q value of r y at root of beam 

S x , S y static mass unbalance associated with x- and 

y-def lections, respectively, mr x , mr y , function 
of z 



x,y 


i, n 


z 

r 


ratio of average static mass unbalance associated 
with x- and y-deflections, respectively, of 
it* 1 interval, to static mass unbalance asso- 
ciated with y-def lection at root 

bending deflections in direction of x- and y-axis, 
respectively, of Cartesian coordinate system 

bending deflection in direction of x- and y-axis, 
respectively, at ith. station 

bending deflection in direction of minor and major 
axes of inertia, respectively, at cross section 
of beam 

distance from root of beam 

(1/ S 2) (I 0 /C 0 ) <ZI x0 /* 0 ) 


length of interval between stations 


coupling coefficient, r JVr ^ 

’ yO ' go 


9 

\ 

c 


angle of torsional vibration 

frequency parameter, latent root of dynamical 
matrices 

dimensionless distance, ^ 
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APPENDIX B 


STATION FUNCTIONS AND DETERMINA.NTAL EQUATIONS 

In the usual influence-coefficient methods for solving dynami- 
cal problems, a continuous body having an infinite number of degrees 
of freedom is replaced by a body having a finite number of degrees 
of freedom. Two principal assumptions are then made, which intro- 
duce inaccuracies into the solution, particularly in the higher 
modes: (l) The inertial loads of all the infinitesimal masses in 

a finite interval can be replaced by a resultant that passes 
through the center of gravity of that interval; and (2) a concen- 
trated load that is the resultant of a distributed lead produces 
the same deflection as the distributed load. 

In order to eliminate these assumptions, Rauscher (refer- 
ence 9) introduced the concept of Station Functions. In this 
approach, it is assumed that the inertial loads and, consequently, 
the deflections are continuous functions along the beam. The 
value of these continuous deflection functions at n reference 
stations must equal the deflections of these reference stations. 

The loading on the beam is therefore a continuous function of the 
deflections of the reference stations. Inasmuch as the deflections 
of the reference stations can be computed from the loading on the 
beam, which in turn is available from the deflections, the deflec- 
tions are therefore obtained as functions of themselves. 

In the method employed here, the deflection of a continuous 
beam is approximated by the sum of a set of continuous Station 
Functions. Each of these functions is chosen: (l) to satisfy the 

boundary conditions of the problem, (2) to vanish at all stations 
other than the one with which it is associated, and (3) to equal 1 
at the station with which it is associated. In terms of these 
functions, the deflections 


e = 


y = 


re 


£ f,(t> 


j=i 
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L,gjU) yj 
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H s,(0 


j=l 


J v V X J 


(Bl) 
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where 

s 


dimensionless distance along "beam, z/5 


e, y, x 

V 7 ) ’ x ) 


torsional deflection, tending deflections in y- 
and x-directions, respectively, functions of £ 

torsional deflection, tending deflections in y- 
and x-directions, respectively, at j'*'* 1 station 


f,( '£), gj(C) 


Station Functions in torsion and tending, respec- 
tively, associated with station 


where 

it (C-J)£ (C-aj.) 

_ Jii 

* (i-j) i (i-a x ) 

J/l 

2 2 (B2) 

« (C-J) £ (£ + a 2^ 4a 3 ) 


it (i-j) i 2 (i 2 +a 2 i-+a ) 

jA 

where rt denotes the product for all values of j except j=i, 

j/i 

and- the constants a-^ ag, and a^ are determined to satisfy the 
boundary conditions at the free end of the team, where the moments 
and the shear forces are zero. 

For a given frequency of vibration, the influence coefficients 
for the system can be derived. The following equations are appli- 
cable at any of the stations: 


fjU) 

gj (£) 
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0, = 5“ | q t ($) / TT ^ 7 d£ 


lo ° Jo c(^) 
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■c 


7± = 5 | q y (0 J fy Ui-tf(i-C) dC a^i + 
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s’ f i 1,7(0 Jo f y (c l‘ ?)<1 ‘ £) 34 4! i ' 


■ 4 ' 1x«) , -xy 


X 
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f„ (£ r O(i-0 at dt 1 - 


s' I 5 X U) I f xy (Sl-SX 1 -?) d ? 4 Si 


— 4 l 


<^(0 t X y ( 5 i^)(H) <34 d^ - 


0 U 0 


6- I V°.L -*y 


f xy (^)(i-0 d{ + 


r 




5 ’ I q x (0 


Ji 


0 UO 

n r»i 


f x (Ci-Od-C) + 


5 <M C) f x (Ci-od-o dCdCj. 


where the loading functions <1^(0, q. x ( £ ) , and q y (£) are 
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q t =' Io?9 + S x <j?x + S y ti) 2 y 

q x = mu£ x + S x O) 2 e (B4) 

2 2 
q y = mti) y + Syi) 0 

where 

C(£^) torsional stiffness of team, function of £ 

I mass moment of inertia per unit length of team, function 

of C 

S x , S y static mass unbalance associated with x- and y-deflec- 
tions, respectively, function of £ 


The integrals of equation (B3) may be written as the sura of 
integrals over each interval by substituting these values for the 
loading functions into equations (B3), using equations (Bl) for 
the deflections, and assuming constant values of m, S x , I, S y , f x , 

f_, and f YV over each interval. The order of summation can then 
be changed to yield 




2 m 0 5 


El 


xO 
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2 “o 5 


El. 


vS Dei J ^ * ‘ r |f yj + cr § 


D 8xlj X 3) 


4 / n 


XO 


2 ^i + 2 D yij yj + 2 D yxij X Jj 
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APPENDIX C 


EXACT SOLUTION 

The differential equations of equilibrium of the cantilever 
beam can be derived from consideration of the forces and the moments 
acting on an elemental section of the beam. These forces and 
moments are shown' in figure 4. 

With sin dqp = dcp and cos dp = 1 and second-order differ- 
entials neglected , from the equilibrium of bending moments 


s 

dz h dz 


V* = 0 


dz 


dcp 

+ Mt — - - V 
* dz ^ 


= 0 


(Cl) 


where 

Ma, lyL bending moment about q- and £-axes, respectively, func- 
^ ^ tion of z 

M torsional moment, function of z 

U 

Va,. V shear force in t- and q~directicns, respectively, func- 
tion of z 

From the equilibrium of forces, • 


dV| 

dz 


£ p 2 

— = mo) S^oj 9 + V^cp’ 


dV. 


^ V p p 

= mu) q + S^co 6 - rp 1 


(C2) 


From the equilibrium of torsional moments, 


2 2 j 2 

- — + Id) 0 + S.oj t + S_uu q = 0 
dz I ^ 


(C3) 
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Upon manipulation of equations (Cl) and (C2) and the substitu- 
tion of the relations 

-"4 


a 2 r, 

*n = EI e^ 


M e = c lf 


> 


(C4) 




where C is the torsional stiffness, a function of z. 

i? 8 ) - (s) ^ i §) - 


a 2 ? 


.2 

a ti 


dz 


=o 


s 

dz 


2 \ EI e 


In 

az 2 y 




■ d 2 q 2 2 2aq> a ( a 2 e > 


a 2 qp a 2 cp 

Ti'T _L_ 


a z ‘ 


EI n 


a z ‘ 


= 0 


a / ae\ 2 2 1 2 

S \ 5/ + I “ e + S s“ £ + V ' 0 

Equations (C5) are the differential equations of equilibrium 
of a cantilever beam vibrating in coupled bending - bending- tors ion 
motion. 


Exact solution. - A uniform y linearly twisted cantilever ‘beam, 
vibrating in coupled bending - bending-torsion motion is considered, 
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Using the change of variables 

Yi * E^ £ 
Y 2 = El, i) 

f 3 

and letting 




Y, = El^r^ 6 




<p t 


‘ *! 


— 7 - 12-36 

El, /ml u. 


-s'- e 

(The value of the constant is taken from reference 11.) 


h 


L e 


£ Vr> 


. ^ 2 /c = /«f _ n = 


m2 /El 


\2 ) 12.2 


4 12 ) 12.36 w 2 - °* 2 ,., 2 £ 1//r t 


0) 


V r 5 s B 


2 , 2 _ . 
r f f r * - € 


£ / r g 

reduces equations (C5) to yield 


Y l' V " VV “ Qy i - 2< P t Y 2 '" “ Qy 3 = 0 


*\ 


Q 


Y 2 ,V - ^tV-^ Y 2 + 2 Vr- RQY 3 = ° 


Y," + Y, + e Q- Y-, + R € -2- Y ? = 0 

3 7 t 3 r t 1 Vb 2 


J 


(C6) 


> (C7) 
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where primes denote successive differentiation with respect to £. 

H M M 

Let Y-j_ = D^e , Y 2 = Ej© } and Yj = Fj© 1 he solutions 
to equation (C7). 


Then 


(A^-cp., 2 Aj 2 - Q) Di - -OFi = 0 

3 4 2 2 

2cpt^i + (A^ -p b A^ - HC1F 1 = 0 


(C8) 


Qd i 


A ' Vt 


^-0 Ei + (A^ + ifi) F t = 0 


p 

A tenth-order equation in X c is obtained by equating the 
determinant of the coefficients of equation (C8) to zero. For the 
case where R = 0, that is, the center of torsion is on . the axis 

of minimum moment of inertia, this equation reduces to a quintic in 


10 2 8 f 4 p 1 

h + + 2< f>t )^i + y>t - q ( i + v^b - 2^ /7 t )j^ 

- U + Vr t )0 "I 4 

+ M + 


Q 


U + + 


7t 


Q c 


cp t 2 (i + l/ 7 b ) <p t 2 e' 
l/ 7 b + - 

L b 7 t 7 t J 


,2 1-c n 3 

+ w =0 

1 Vb 


(C9) 


Let the roots of equation (C9) be A^, A ^ A^, A g, A^, A 


A„ , A . , A _ , and A , where 
4* -4 * 5* -5 7 


h = - A_^, i = 1, 2, . . .5 


Then 
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where 


Y l = 


Y 2 = 


5 V 

t V 

i=-5 

i /0 




t *!■ 




> (CIO) 


i=-5 

i/0 


Y 3 = 


A , K ±l 

2_, F i e 
i=-5 

i^O 


J 


E* = 


• 2c Pt^i 




1 " (x/.^-n/r!,) 


Dj H GjDi 


Fl ‘ Di HlDl 

r t \ +n 

i = ±1, ±2, ±3, ±4, ±5 


(Cll) 


The boundary conditions are 
at £ = 0 

Y (o) = Y '(0) = Y (0) = Y *( 0 ) = Y (0) = 0 
at-C -1 (C12) 

Y^'Cl) = Y 1 "' (1) = Y 2 "(l) = Yg'" (1) - Y 3 '(l) - 0 . 


The substitution of equation (CIO) into the appropriate equa- 
tions of (C12) will yield 10 homogeneous equations in the 10 vari- 
ables Di. For a nontrivial solution, the determinant of the coef- 
ficients must vanish. Thus: 
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in 

in ,< in 
O in W 
o 


to 

to l 


o’ * 



rH <<? 


© 

CM 

rH 


H 

<< ^ "h W ri 00° 

© 

to 

to 


& H 

,< 

i-H 

,< 

& 

o 

c. 


(C13) 
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The process of solution is as follows: The value of the 

determinant in equation (C13) must he -plotted against frequency; 
the value of frequency for which this determinant becomes zero is 
thereby obtained. This procedure involves first solving the quin- 
tic equation ( C 9 ) for each assumed value of frequency parameter Q. 
and then calculating the elements of the determinant from equa- 
tion (Oil). The value of frequency for which the determinant of 
equation (C13) equals zero is the natural frequency of vibration 
of the beam. 

The exact solution is obviously long and laborious. Numerical 
calculations have been made for the case where there is no tor- 
sional coupling with the bending modes. In this case, equation (C12) 
may be replaced by a biquadratic equation in A-j , and the order 
of the determinant of equation (C16) is reduced to 8. The results 
of these calculations are included in the body of the report, where 
a comparison is made with values given by the approximate solution. 
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APPENDIX D 


ELEMENTARY PROPERTIES OF MATRICES 

A matrix is a set of numbers or other elements arranged in 
rows and columns that obey certain rules of addition and multipli- 
cation. Thus 


G 

is a matrix of numbers containing two rows and three columns. 

" 4 “ 

1 

0 

_- 2 _ 

is a matrix containing four rows and one column. It is called a 
column matrix. Similarly, a set of numbers arranged in a single 
row is designated a row matrix. A nXm matrix is one containing 
n rows and m columns. 

Let an nXm matrix be designated A. Then the element in 
the i bb row and column of A will be designated a^j 

(the first subscript designating the row, the second subscript the 
column) . 

The operations of addition and multiplication of matrices are 
now defined as follows: 

If a matrix A • is added to a matrix B (B must have the 
same number of rows and columns as A), the result is a matrix C 
given by 


0 3 

-4 9 


] 


C = A + B where c^j = a^j + 

If an nXm matrix A is multiplied by a matrix B, the 
result is a matrix C given by 


C = AB where c^j 



a ik b kj 



30 


NACA TN 2300 


that is, the element in the i^’ n row and j ua column of C is 
obtained by taking the sum of the products obtained by multiplying 
the elements of the i^ row of A with the corresponding ele- 
ments of the column of B. 


<th 


It is to be noted that the order of multiplication is impor- 
tant, that is, in general, BA / AB, and that in order to multiply 

two matrices together, the number of columns of the first must 
equal the number of rows of the second. 

Examples : 

D °J-[i -3 ■ [i 3 
i 3 [-; .3 ■ ;] 

Pi 7~| [l o| [l5 35l 

15 -2J [2 5j = 1-7 -lOj 

For a complete discussion of the properties of matrices, see 
reference 12 . 


REFERENCES 

1. Myklestad, N. 0.: Vibration Analysis. McGraw-Hill Book Co., 

Inc., 1944. 

2. Targoff, Walter P.: The Associated Matrices of Bending and 

Coupled Bending-Torsion Vibrations. Jour. Aero. Sci., vol. 14, 
no. 10, Oct. 1947, pp. 579-582. 

3. Fettis, Henry E.: The Calculation of Coupled Modes of Vibration 

by the Stodola Method. Jour. Aero. Sci., vcl. 16, no. 5, 

May 1949, pp. 259-271. 

4. Mendelson, Alexander, and Gendler, Selwyn: Analytical Deter- 

mination of Coupled Bending-Torsion Vibrations of Cantilever 
Beams by Means of Station Functions. NACA TN 2185, 1950. 



NACA TN 2300 


31 


5. Morris, J., and Head, J. W. : Coupled Flexural Vibrations of 

the Blades of a Propeller and Torsional Vibration of an Engine 
Crankshaft System. R. & M. No. 2011, British A.R.C., Oct. 
1942. 

6. Love, E. R., and Silberstein, J. P. 0.,. assisted by Radok, 

J. R. M.: Vibration of Stationary and Rotating Propellers. 

Rep. SM. 48, Div. Aero., Council Sci. Ind. Res. (Australia), 
May 1946. 

7. White, Walter T.: An Integral-Equation Approach to Problems of 

Vibrating Beams. I. Jour. Franklin Inst., vol. 245, no. 1, 
Jan. 1948, pp. 25-38; II, vol. 245, no. 2, Feb. 1948, 
pp. 117-133. 

8. Plunkett, R.: A Matrix Method of Calculating Propeller-Blade 

Moments and Deflections. Jour. Appl. Mech. , vol. 16, no. 4, 
Dec. 1949, pp. 361-369. 

9. Rauscher, Manfred: Station Functions and Air Density Varia- 

tions in Flutter Analysis. Jour. Aero. Sci., vol. 16, no. 6, 
June 1949, pp. 345-353. 

10. Geiger, J.: Influence of Twist in Turbine Blades on Natural 

Frequencies and Directions of Vibrations. Eng. Digest, 
vol. XI, no. 4, April 1950, pp. 115-118. 

11. Timoshenko, S.: Vibration Problems in Engineering. D. Van 

Nostrand Co., Inc., 1928, p. 344. 

12. Frazer, R. A., Duncan, W. J., and Collar, A. R.: Elementary 

Matrices. Cambridge Univ. Press (London), 1938. 



32 


NACA TN 2300 


TABLE I - STATION NUMBERS, n = 1 


M Jk N Jk p Jk 


Q 




^5/isJ [2/3J [5/12] [3/2 oj |j/ieoj 

L 'jk M ’jk N ’jk , p ’jk Q'.1k 

[i.3/4 5^ ^2/5^ |l3/ 45J |7l/630j j3l/l008j 


TABLE II - STATION NUMBERS, n = 2 


M 


N 


8/15 -31/240"! pLl/12 -13/48] p/15 -31/24(5 

7/60 94/24oJ L 5 / 12 29/48 J |_ 8 / 15 239/24Q 


[ 


P Q 

.183333 -0.037500"] fo. 046032 -0. 008135"] 

.025000 .143750J L. 029365 .181448J 


L’ 


t 


.224675 .3100111 I .627273 


JL 


M' 


.367100 -0.0348751 IS. 5363 64 -0.060795' 


5 1 


,448674| 


N’ 




0.367100 -0.0348751 fO. 137933 -0.011252" 

_ .851948 .758685J L .057955 .118463 


Q’ 

[ 0.036616 -0.002614*1 

. .069733 . 150415 J 
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TABLE III - STATION NUMBERS, n = 3 
L M 


0.545833 

-0.241667 

0.106019 


0. 950000 

-0.525000 

0.235185 

. 137500 

.450000 

-.078240 


. 450000 

.725000 

-. 153704 

-.020833 

.141667 

.376388 


-. 050000 

.475000 

.568519 


N 


1 


P 


0. 545833 

-0.241667 

0.106019 


0. 186310 

-0.068452 

0.029563 

.587500 

1.175000 

-. 231944 


. 032143 

.160714 

-. 023677 

-.120833 

1.091667 

1.513426 


-. 005357 

.031548 

. 139749 


r* 


Q. 


“i 



0.046577 -0.014583 0.006222 
.038244 .202083 -.028963 

-.011756 .068750 .316408 



L' 



M' 


0.399646 

-0.083406 

0.022325 


0.596268 

-0. 149356 

0.040630 

. 175000 

.391964 

-. 038816 


. 533205 

.602896 

-. 072928 

-. 045142 

. 224317 

.308509 


-. 097426 

.625418 

.445812 


N' 



P' 


0.399646 

-0.083406 

0.022325 


0.148013 

-0. 026378 

0.006972 

.708205 

.994660 

-.111744 


.042560 

. 143948 

-.012081 

-. 239994 

1.475153 

1.200133 


-.012798 

. 057937 

.116007 


r . 


Q' 


“1 

'\NACA^' 


0.038884 -0.006034 0.001579 
.050843 .181698 -.014830 

-.028318 .127659 .267865! 
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TABLE IV - STATION NUMBERS, n = 4 


L M 


0.576455 
. 128307 
-.023545 
p .009789 

-0. 292657 
.459226 
.169643 
-.036606 

0.229365 
-. 126190 
.407144 
.162699 

-0. 096544 
.045874 
-. 055919 
,365434 

' 

1.022222 
.429630 
-.051852 
_. 022222 

-0.647917 

.747917 

.518750 

-.085417 

0.522822 
-.255556 
. 633333 
. 522222 

-0.221701 

.094650 

-.105961 

.543924 

0.576455 

.557937 

127249 

.076455 

N 

-0.292857 

1.207143 

1.207143 

-.292857 

0.229365 

-.381746 

1.673810 

1.729365 

-0.096544 

.140724 

-.267841 

1.997206 


0. 194478 
.029597 
~. 006581 
.002612 

P 

-0.081920 

.163021 

.041295 

009598 

0. 06 2798 
-. 037401 
.148512 
.037202 

-0.026267 
.013391 
— • 017322 
. 136803 


Q, 


0.048240 
.035167 
-.014547 
. 008359 


-0.017295 
. 204826 
. 090642 
-. 030688 


0.013040 
-. 045624 
.335791 
.118397 


-0.005428 

.016301 

038574 

.446729 


L* 


M 1 


0.413738 
. 164798 
-. 037937 
.020126 


-0. 116662 
.424050 
. 181143 
-. 052641 


0.065879 

086085 

.381486 

.233289 


-0.019023 
.021577 
-. 034723 
.305180 


0.623188 
.511682 
-. 082891 
.042276 


-0.211987 
.667412 
. 544025 
-.112648 


0.122052 
-. 166738 
• 582158 
.643846 


-0.035456 
• 042626 
-.064723 
.438962 


N» 


P» 


0.413738 
. 676680 
203719 
. 146954 


-0. 116662 
1.091462 
1.269193 
390585 


0. 065679 
252823 
1.545802 
2.164827 


— 0. 019023 
. 064205 
-.164169 
1.622066 


0.152256 

.039616 

-.010651 

.005795 


-0.036502 
. 153469 
.044508 
-.015000 


0.020304 

026235 

.140822 

.060554 


-0.005836 
. 006481 
-.010869 
.117037 


Q* 



0.039818 
.047267 
-.023551 
. 018630 


-0.008281 
.193310 
.097745 
-. 048203 


0.004551 
-. 032114 
.318707 
.194016 


-0.001303 

.007917 

024222 

.382752 



TABLE V - STATION NUMBERS 
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Figure 2. — Uniform naturally twisted cantilever beam. 
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